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I.  INTRODUCTION 


A typical  problem  in  the  measurem«>nt  of  collision  cross  sections  is 
the  determination  of  resonance  peaks  of  narrow  width.  Often,  the  energy  re- 
solving power  of  the  apparatus  is  inadequate  to  define  clearly  the  shape  of  the 
resonance  peak.  When  an  electron  or  molecular  beam  is  used,  the  quantity 
that  is  measured  is  an  energy-averaged  cross  section  a that  is  related  to  the 
true  cross  section  a by 


a(V^)  = 


Va(V)f(V;  V 


t)  dV 


t)  dV 


(1) 


where  V is  the  energy  of  a particle  in  the  beam,  is  a characteristic  energy 
of  the  beam,  t is  a characteristic  energy  spread  of  the  beam,  and  f(V;V^,  t) 
is  the  energy  distribution  of  the  beam.  If  the  energy  spread  of  the  beam  is 
negligible  compared  to  V^, 

The  measurement  of  the  distribution  function  f is  generally  very  dif- 
ficult. For  this  reason,  the  approximation  that  o = o sometimes  is  made. 

If  the  energy  spread  of  the  beam  is  made  sufficiently  small,  this  approxima- 
tion obviously  will  be  tolerable.  However,  in  many  cases  of  practical  interest 
(Refs.  1 and  2),  the  energy  spread  cannot  be  reduced  sufficiently  to  justify 
the  use  of  this  approximation.  In  these  cases,  even  an  approximate  estimate 


^P.  J.  Chantry,  Attachment  Measurements  in  Halogen  Bearing  Molecules, 
Westinghouse  R&D  Center  (March  1978). 

^P.  Mahadevan,  M.  Fpstein,  and  R.  Holland,  Jr.,  Final  Report;  Plasma 
Cross-Section  Measurements  Relevant  to  Atmospheric- Pressure  High-Pqwe^ 
DF/HF  Chain  Lasers.  ATR- 78(7585)- 1,  The  Aerospace  ::orporation, 

El  Segundo,  California  (31  August  1978). 
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of  the  distribution  fiinction  could  result  in  a significantly  better  estimate  of 
the  true  cross  section  if  the  integral  equation,  Eq.  (1),  can  be  solved. 

In  principle,  the  solution  of  Eq.  (1)  is  relatively  straightforward. 

Several  numerical  techniques  have  been  devised,  most  of  which  are  variations 
of  a collocation  technique  (Ref.  3).  In  this  approach,  it  is  assumed  that  ex- 
perimental measurements  have  been  made  at,  e.g.,  N values  of  V . N values 

di 

of  V at  which  a is  desired  are  then  chosen.  The  function  a is  interpolated 
analytically  between  adjacent  vcilues  of  V so  that  by  writing  Eq.  (1)  N times 
(at  the  N values  of  V ),  a set  of  N algebraic  equations  are  obtained  for  the  N 
unknown  values  of  a. 

Unfortunately,  the  numericad  computations  associated  with  this  approach 
are  difficult  to  perform  because  the  matrices  that  have  to  be  inverted  are 
ill-conditioned.  The  accuracy  of  the  method  is,  therefore,  rather  poor  unless 
the  data  are  very  smooth.  This  condition  is  not  likely  to  be  satisfied  for  the 
case  of  a narrow  resonance  peak. 

Some  researchers  have  recognized  this  problem  and  have  approached 

the  problem  of  smoothing  the  data  with  a very  sophisticated  technique  (Ref.  4). 

In  many  of  the  cases  of  interest,  the  use  of  this  type  of  technique,  provides 

significantly  improved  results.  However,  since  this  sophisticated  approach 

involves  the  high-frequency  filtering  of  the  data  (to  eliminate  experimental 

scatter),  it  may  still  fail  to  produce  adequate  results  when  the  width  of  the 

resonance  peak  is  as  small  as  a few  times  the  uncertainty  in  the  resolution 

of  V . 
a 

A numerical  approach  is  described  here,  which  appears  to  be  useful 
when  the  cross  section  is  rapidly  varying  over  some  energy  range,  but  slowly 
varying  over  an  energy  range  large  compared  to  the  energy  spread  of  the  beam. 
An  exact  solution  of  the  integral  equation  is  also  given  for  a form  of  the  distri- 
bution function  that  is  of  practical  interest. 

_ 

K.  Berkling,  R.  Helbing,  K.  Kramer,  H.  Pauly,  C.  H.  Schlier,  and 
P,  Toschek,  Z.  Phys.  166,  406  (1962). 

4 

J.  D.  Morrison,  J.  Chem.  Phys.,  39,  ZOO  (1963) 
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n.  ANALYSIS 


It  is  assumed  that  the  true  cross  section  has  the  qualitative  form  shown 
in  Fig.  1,  i.  e.  , a narrow  resonance  peak  at  V = 0 and  smooth  behavior  at 
large  values  of  V.  The  method,  in  fact,  is  applicable  to  cross-section  fvmc- 
tions  containing  multiple  peaks  of  arbitrary  locations  provided  the  measure- 
ments extend  over  an  energy  range  that  is  large  compared  to  the  sum  of  the 
widths  of  the  resonance  peaks  and  the  energy  spread  t of  the  beam.  This 
method  has  been  developed  specifically  for  two  particular  forms  of  the  beam- 
energy  distribution.  These  results  are  sufficient  to  indicate  how  the  method 
can  be  applied  to  other  energy  distributions. 

The  particular  beam-energy  distribution  considered  is  appropriate  for 
many  experimental  investigations  of  cross  sections  for  electron  heavy-particle 
collisions.  If  the  t>eam  is  formed  by  emitting  electrons  from  a cathode  and 
accelerating  them  through  an  anode  aperture  to  form  an  electron  beam,  the 
distribution  function  can  be  reasonably  approximated  by 

0 V < Va 

(V  - V )/T 

® V > V 

a 


where  V is  the  accelerating  potential  of  the  anode  with  respect  to  the  cathode, 

a 

and  T is  the  energy  spread  of  the  beam.  For  this  form  of  f,  Eq.  (1)  becomes 


7(V^)  = 


t(t+  V. 


0 


Va(V)  e 


-(V 


V )/T 
a 


dV 


(3) 


Before  the  numerical  technique  that  has  been  developed  is  described, 
an  exact  solution  of  Eq.  (3)  for  the  special  case  of  constant  t is  discussed 
briefly.  The  quantity  t is  a constant  when  the  energy  spread  of  the  beam  is 
independent  of  the  accelerating  potential.  The  exact  solution  is  obtained 
almost  trivially  by  differentiating  both  sides  of  Eq.  (3)  with  respect  to  V and 
solving  the  result  for  a.  The  solution  for  constant  t is 


a(VJ  =a{VJ  - 


t(t  + V ) 

^ 

V 

a 


do 

dV 

a 


(4) 


In  practice,  this  method  has  certain  distinct  advantages  and  disadvantages 
The  primary  advantage  is  that  no  elaborate  computations  are  needed  to  calcu- 
late a (if,  of  course,  t is  assumed  to  be  constant  and  adequately  known).  On 
the  other  hand,  the  derivative  ofo  must  be  calculated,  and,  consequently, 
scatter  in  the  data  could  seriously  affect  the  accuracy  of  the  results.  The 
problem  of  scatter  can  be  reduced  somewhat  by  one  of  two  techniques.  In  one 
technique,  the  data  is  smoothed.  Relatively  unsophisticated  smoothing  algo- 
rithms should  be  adequate  to  provide  a substantial  improvement  in  the  esti- 
mate of  a.  In  the  other  technique,  the  scatter  often  becomes  large  when  the 
experimental  apparatus  is  pushed  to  the  limits  of  its  resolving  capability.  A 
beam  with  a larger  energy  spread  tends  to  be  easier  to  deal  with  and  has  cor- 
respondingly less  scatter.  Therefore,  if  the  parameter  t can  be  defined  with 
sufficient  accuracy,  a beam  with  an  extremely  small  energy  spread  is  not 
needed. 

Note  that  near  the  peak  of  the  measured  cross  section,  Eq.  (4)  is  in- 
determinate. By  the  application  of  L'Hospital's  rule. 


is  obtained.  The  difficulty  of  computing  acceptably  accurate  second  derivatives 
from  experimental  data  would  appear  to  limit  the  usefulness  of  this  approach 
in  obtaining  the  maudmum  value  of  the  resonance  peak.  However,  for  moderate 
values  of  T,  it  may  be  possible  to  minimize  scatter  so  that  the  curve  of  o'  ver- 
sus is  relatively  smooth,  and  a tolerable  estimate  of  d^a/dV^  at  the  origin 
can  be  made.  Various  types  of  data  smoothing  techniques  have  been  developed 
that  provide  estimates  for  upper  and  lower  bounds  for  the  derivatives  of  the 
function  being  fitted,  which  can  be  used  to  estimate  upper  and  lower  bounds 
for  the  peak  value  of  a. 

For  the  more  general  case  of  nonconstant  t,  it  is  necessary  to  resort 
to  a numerical  technique  for  the  solution  of  Eq.  (3).  The  present  technique 
requires  that  the  energy  range  covered  by  the  experiment  be  large  compared 
to  the  sum  of  the  widths  of  the  resonance  peaks  appearing  in  the  energy  inter- 
val of  interest,  2uid  large  with  regard  to  t . In  this  case,  there  will  be  an 
energy  range  in  which  a because,  in  that  range,  the  distribution  function 
is,  to  first  order,  a 6 function  compared  to  the  cross  section  curve.  It  is 
most  convenient  if  that  region  occurs  at  energies  that  are  large  compared  to 
the  energy  at  which  the  resonance  peak  occurs. 

For  such  conditions,  the  value  of  a at  large  energies  is  measured 
directly,  i.  e.,  a = a for  greater  than  same  value,  e.g.,  VaQ*  For  a value 
of  V^,  slightly  less  than  V^q,  e.g.,  and  if  sufficiently  small, 

a can  be  approximated  in  the  interval  Vai  ^ V < Vao  by  means  of 

•a(V)=ao+^J-r-^  (V-  Vq)  (6) 
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t' 
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where  the  approximation  a=  a for  V > Vg^Q  has  been  used  in  the  second  integral 
in  Eq.  (7).  If  the  interest  is  centered  on  determining  the  magnitude  and  shape 
of  the  resonance  peak(s),  the  second  integral  will  be  small,  and  a negligible 
error  will  be  incurred  if  ^ is  approximated  in  any  suitable  way  for  V greater 
than  the  energy  range  covered  in  the  experiment.  If  is  an  energy  at  which 
a has  been  measured,  the  only  unknown  in  Eq.  (7)  is  s o(Va^j).  The  first 
integral  in  Eq.  (7)  can  be  evaluated  analytically,  and  the  result  is  a linear 
algebraic  equation  for  a^. 

When  Oj  is  determined,  the  same  process  can  be  repeated  at 
a value  slightly  less  than  In  this  way,  a marching  type  of  calculation 

can  be  constructed,  which  permits  the  calculation  of  at  each  all  the 
way  to  V = 0,  with  one  linear  algebraic  equation  solved  at  each  step.  The 
result  for  a.  = a(Vj)  is. 


where 


k = 


0 i = 1 

1 i = 2,  3,  •••N 


A..  = V.  .C..  - D.. 
ij  J-1  iJ  iJ 


B..  = V.C..  - D.. 
ij  3 IJ  13 


(9) 

(10) 

(11) 


-(V,/t  ) .(V  Jr  ) 

V.  , - V. 

3-1  3 


(12) 


2 2 2 2 -^^i-l/V 

(2t‘ + 2t.V.  + Vp  e J -(2Tf + 2TjV._j  + Vpj)e  ^ " 

^ij  ^ V.  , - V. 

•'  3-1  3 


(13) 


and  N + 1 is  the  number  of  values  of  V at  which  a is  measured. 


At  the  i'th  step,  a can  be  calculated  directly  from  the  value  of  7 at  the 
i'th  step  and  from  the  values  of  a already  calculated  in  the  previous  steps. 
The  calculation,  therefore,  is  extremely  rapid  and,  as  shown  in  the  follow- 
ing section,  accurate. 
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m.  ILLUSTRATIVE  EXAMPLE 


The  usual  approach  to  validating  the  technique  developed  here  would  be 
to  perform  the  calculation  for  a set  of  experimental  data  and  to  compare  the 
results  with  the  results  of  calculations  obtained  with  other  methods.  This 
approach  has,  however,  several  drawbacks.  First,  no  suitable  data  are 
available  in  this  laboratory.  Second,  of  the  data  available  in  the  literature, 
the  experimental  conditions  are  generally  not  given  in  sufficient  detail  to  per- 
mit an  unambiguous  comparison  to  be  made.  Further,  in  order  to  compare 
with  other  calculation  methods,  those  methods  must  be  available.  The  more 
sophisticated  (and  therefore  presumably  more  accurate)  methods  involve 
moderate-to-long  computer  codes.  Even  if  an  effort  to  develop  those  codes 
were  made,  the  results  of  the  comparisons  would  not  be  conclusive  since 
none  of  the  methods  are  exact.  Comparisons  with  published  calculations  are 
generally  inconclusive  because  details  of  the  smoothing  are  usually  not 
available . 

Because  of  these  drawbacks,  another  approach  was  taken.  It  would  be 
ideal  to  apply  the  calculation  technique  to  experimental  data  to  estimate  the 
true  cross  section  and  then  compare  the  result  with  the  real  true  cross  sec- 
tion. Since  the  real  true  cross  section  is  never  known  exactly,  the  calcula- 
tion was  carried  out  as  follows:  An  analytically  simple,  but  qualitatively, 
realistic  form  for  the  true  cross  section  was  assumed  and  used  in  a thought 
experiment.  It  was  assumed  that  an  electron  beam  with  a given  energy  dis- 
tribution interacted  with  the  gas  in  question,  and  measurements  were  made 
of  the  energy  averaged  cross  section.  If  the  measurements  were  made  with 
infinite  precision,  the  measured  cross  section  could  be  obtained  by  quadrature 
from  Eq.  (1).  A perturbation  of  the  exact  value  of  a was  generated  at  each 
value  of  by  Monte  Carlo  techniques  in  order  to  simulate  the  unavoidable 
experimental  scatter.  These  values,  which  include  scatter  with  a specified 
vjtlue  of  standard  deviation,  were  then  treated  as  raw  experimental  data.  The 


comparison  of  the  calculated  true  cross  section  with  the  true  cross  section 
from  which  the  "experimental"  data  were  derived  provides  a direct  evaluation 
of  the  accuracy  of  the  method. 

Accordingly,  the  form  for  the  true  cross  section  that  was  selected  was 


o =aoe 


-(V/V^) 


where  is  the  peak  value  of  the  cross  section  and  is  a measure  of  the 
width  of  the  resonance  peak.  By  defining 


V 

V = 

a 


S = -2- 


and  by  means  of  Eq.  (1)  with  the  distribution  function  given  by  Eq.  (2), 


t + (1  + t)v 


(1  +t)‘'(t  + v) 

is  obtained.  Note  that  this  result  is  true  for  die  general  case  t = t(v). 

The  thought  experiment  is  then  performed  under  the  assumption  that, 
at  each  value  of  v,  the  "measurement"  is  sampled  from  a normal  distribution 
with  a meem  of  S(v)  and  a standard  deviation  that  is  a fixed  fraction  A of  ^(v). 

Calculations  were  made  for  t = 1.  Thirty  three  values  of  v were  selected 
between  v = 0 and  v = 10.  At  each  value  of  v,  three  "measurements"  were 
made  for  a total  of  99  data  points.  Approximately  random  numbers  were  ob- 
tained by  selecting  the  last  four  digits  in  a tabulation  of  natural  logarithms 
given  to  five  places  (Ref.  5).  The  first  of  those  four  digits  was  used  to 


*R.  S.  Burington,  Handbook  of  Mathematical  Tables  and  Formulas.  Handbook 
Publishers,  Inc.,  iSandusky,  Ohio  U953)  p.  2^1. 


determine  the  sign  of  the  deviation  from  the  mean  (plus  for  an  even  digit  and 
minus  for  an  odd  digit).  Any  minor  deviation  from  a strictly  random  selection 
of  numbers  can  be  interpreted  as  modeling  the  bias  in  am  experiment  and  thus 
adds  to  the  analogy  to  a real  experiment. 

The  results  of  numerically  simulating  an  experiment  in  this  manner  are 
shown  in  Fig.  2 for  6 = 0. 1.  The  true  cross  section  and  the  exact  energy- 
averaged  cross  section  are  shown  for  comparison.  For  this  particular  reali- 
zation of  the  "e3q>erimental"  data,  all  of  the  data  points  are  within  ±20%  of 
the  exact  energy- averaged  cross  section. 

In  a real  experiment,  it  would  be  appropriate  to  smooth  the  data  if  it 
is  believed  that  there  is  no  fine  structure  that  might  be  lost.  The  smoothing 
interval  is  generally  chosen  with  this  in  mind.  Accordingly,  the  simulated 
data  were  subjected  to  a simple  smoothing  routine,  by  fitting,  in  a least- 
squares  sense,  successive  groups  of  15  points  by  parabolas.  However,  this 
crude  fitting  technique,  constitutes  a two-edged  sword.  On  one  edge,  the 
scatter  is  eliminated;  on  the  other,  in  regions  where  the  scatter  is  small, 
an  error  is  introduced  because  the  true  cross  section  cannot  be  fit  exactly  by 
a parabola. 

The  exact  energy- averaged  cross  section  was  subjected  to  the  same 
smoothing  procedure,  and  the  resultant  data  were  used  to  calculate  a with  the 
use  of  Eq.  (8)  in  order  to  assess  the  sensitivity  of  the  calculation  technique 
to  errors  induced  by  smoothing.  The  results  of  diis  calculation  are  shown 
in  Fig.  3.  The  error  introduced  by  smootiiing  is  negligible,  except  very 
close  to  V = 0.  Careful  consideration  of  Eqs.  (8)  through  (13)  indicates  that 
this  is  because  of  the  finite  difference  nature  of  the  calculation  near  v = 0 
rather  dian  the  existence  of  a peak  in  the  cross  section.  From  these  results, 
it  appears  that  the  calculation  of  resonance  peaks  located  at  least  20%  of  the 
energy  spread  of  the  beam  away  from  zero  would  not  be  affected  by  the  smooth- 
ing process. 


The  true  cross  section  calculated  from  the  data  of  the  simulated 
exp>eriment  is  shown  in  Fig.  4.  There  is  good  agreement  except  for  v less 
than  about  0.  5.  A more  sophisticated  smoothing  technique  would  probably 
have  significantly  improved  the  results  for  small  v.  This  is  borne  out  not 
only  by  the  results  but  also  because  the  exact  solution,  the  dashed  curve  from 
Eq.  (4),  also  indicates  that  there  are  marked  irregularities  for  small  v, 
indicating  that  the  function  is  being  inadequately  fit  in  that  region. 

As  another  example  of  this  same  general  approach,  the  new  method  was 
applied  to  Eq.  (1)  for  another  special  distribution  function  of  importance.  One 
of  the  more  common  techniques  for  making  measurements  of  the  type  under 
consideration  here  is  the  retarding  potential  difference  (RPD)  method  (Ref.  6). 
In  this  method,  a narrow  slice  of  the  energy  distribution  of  the  beam  is  selected 
in  an  effort  to  increase  the  monochromicity  of  the  beam.  As  a first  approxi- 
mation, it  is  often  assumed  that  the  distribution  function  is  represented  by 


f = 


constant 

0 


V < V < V + € 
a a 

elsewhere 


(17) 


By  means  of  this  expression  in  Eq.  (1), 

* V, 

a 

is  obtained. 

The  calculation  proceeds  as  previously  described,  i.  e.  , it  is  begun  in  a 
region  where  a "■  a.  For  example,  for  the  situation  corresponding  to  Fig.  1, 
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the  calculation  is  begun  at  large  V^.  Then,  by  taking  a small  step,  AV  to 
- AV,  and  approximating  o(V)  by  a linear  function  in  that  interval, 

Eq.  (18)  can  be  integrated  to  provide  an  explicit  expression  for  o(V  - AV) 

cl 


a(V^  - AV)  = 


AV(V^  - AV  + "I  «)7(V^  - AV)  - « -I  - AV)  + € 


AV(V^  - AV)  -€(-|  V^  - Av) 


a(V  ) 

St 


(19) 


At  V^  - AV,  a is  given  in  terms  of  the  energy-averaged  cross  section  at  that 
point,  the  true  cross  section  at  V (which  is  known),  and  ihe  width  of  the 
energy  "slice"  c.  By  repeating  this  process  and  marching  inward  to  lower 
values  of  V , the  complete  function  a(V)  can  be  determined. 

Equation  (19)  has  been  applied  to  the  recent  data  acquired  by  Chantry 
for  the  dissociative  attachment  of  electrons  to  F,  (Ref.  1).  At  small  energies, 
the  true  cross  section  is  observed  to  be  about  a factor  of  two  greater  than  the 
measured  cross  section.  In  performing  these  calculations,  the  tabulated  data 
presented  by  Chantry  were  used  without  smoothing  (but  modified  at  low  ener- 
gies as  discussed  by  Chantry).  The  apparent  maximum  at  nonzero  energy  in 
Fig.  5 may  be  an  artifact  related  to  (he  limited  resolution  of  the  apparatus  or 
to  the  sensitivity  of  the  calculation  technique  to  small  errors  near  V = 0,  or 
both.  Regardless,  the  smoothness  of  the  results,  which  demonstrates  the 
stability  of  the  calculation,  is  apparent.  Note  that  this  calculation  was  easily 
performed  with  a hand  calculator. 
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IV.  CONCLUSIONS 


A relatively  simple  calculational  technique  has  been  devised  to  solve 
the  integral  equation  that  relates  the  experimentally  measured  cross  section 
to  the  true  cross  section.  This  technique  is  relatively  free  of  the  problems 
encountered  in  previous  methods,  viz,  the  stability  of  the  calculation  technique 
and  consequent  sensitivity  to  scatter  in  the  data.  Scatter,  of  course,  did 
affect  the  present  results,  but  moderate  amounts  of  scatter  do  not  destroy 
the  stability  of  the  calculation.  In  fact,  in  the  numerically  simulated  experi- 
ments, if  the  solution  is  smoothed  by  passing  a curve  through  the  midpoints 
of  the  linear  segments  constituting  the  solution,  a nearly  exact  solution  is 
obtained.  This  smoothing  of  the  solution  was  applied  in  the  second  example 
described  in  Section  III.  Smoothing  of  the  raw  data  generally  results  in  a 
smoother  solution,  just  as  with  any  of  the  other  techniques. 

If  the  energy  distribution  of  the  beam  is  poorly  known,  this  uncertainty 
will  carry  over  into  the  solution.  On  the  other  hand,  unless  the  beam  is  nearly 
monoenergetic  compared  to  the  width  of  the  resonance  peak  of  the  cross  sec- 
tion, a significant  improvement  in  the  estimation  of  the  true  cross  section  can 
be  expected  by  applying  the  present  deconvolution  technique. 

For  the  important  and  practical  case  of  a translated  Maxwellian  beam 
energy  distribution,  an  exact  solution  has  been  found  for  the  case  of  constant- 
beam  temperature.  Since  the  result  contains  the  derivative  of  the  measured 
cross  section,  the  numerical  accuracy  of  results  obtained  from  this  solution 
are  limited. 
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U^BORATORY  OPERATIONS 

Th«  Laboratory  Operationa  of  The  Aeroapaca  Corporation  la  conducting 
exparlmantal  and  theoretical  inveatigationa  neceaaary  for  the  evaluation  and 
application  of  acientiRc  advancea  to  new  military  concepta  and  ayatema.  Var- 
aatility  and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
peraonnal  in  dealing  with  the  many  problema  encountered  in  the  nation'a  rapidly 
developing  apace  and  miaaile  ayatema.  Expertiae  in  the  latent  acientiflc  deval- 
opmenta  ta  vital  to  the  accompliahment  of  taaka  related  to  theae  problema.  The 
laboratoriaa  that  contribute  to  thia  reaearch  are: 

Aero£hyatcaL«boratorj[;  Launch  and  reentry  aerodynamlca,  heat  trana- 
fer,  reenti^  phyaica.  chemical  kinetica,  atructural  mechanica,  flight  dynamical 
atmoapherlc  pollution,  and  high-power  gaa  laaera. 

Chemiatrv  and  Phvalca  Laboratory:  Atmoapharic  reactiona  and  atmon- 
pherlc  optica,  chemical  reactiona  in  poHuted  atmoapharea,  chemical  reactiona 
of  axcit^  apeciaa  in  rocket  plumea,  chemical  thermodynamic  a,  plaama  and 
laaer-lnducod  reactiona,  laaer  chemtatry,  propulaion  chemiatry,  apace  vacuum 
and  radiation  effecta  on  materiala,  lubrication  and  aurface  phenomena,  photo- 
aenaitiva  materiala  and  aennora,  high  preciaion  laaer  ranging,  and  the  appli- 
catloc  of  phyaica  and  chemiatry  to  problema  of  law  enforcement  and  biomedicine. 

Eleetronlca  Reaearch  Laboratory:  Electromagnetic  theory,  devlcea,  and 
propagation  pkenomana,  including  pkaarna  electromagnetica;  quantum  electronica, 
laaera,  and  electro-optica;  communication  aciancea,  appliad  alactronica,  aami- 
conductlng,  auparconducting,  and  cryatal  davica  phyaica,  optical  and  acouatical 
Imaging;  atmoapharic  pollution;  millimeter  wave  and  far-infrared  technology, 

Materiala  Sctencea  Laboratory;  Development  of  new  materiala;  motel 
matrix  compoaitea  and  new  forma  of  carbon;  teat  and  evaluation  of  grapUta 
and  ceramlca  in  reentry;  apacecraft  materiala  and  alactronic  componanta  in 
nuclear  weapona  environment:  application  of  fracture  mechanica  to  atreaa  cor- 
roalon  and  fatigua -induced  fracturee  in  atructural  matala. 


e 


I . Space  Sctencea  l4>boratery;  Atmoapherlc  and  ionoapharic  phyaica,  radia- 

f tlon  tram  the  atmoapharo,  danalty  and  compoaition  of  flia  atmoai^ra,  aurorae 

I and  airglow;  magnatoapheric  phyaica,  coamic  raya,  ganeratlon  and  propagation 

: of  plaama  wavea  in  the  magnetoaphera;  aolar  phyaica,  atudiaa  of  aolar  magnetic 

j fielda;  apace  aatronomy,  x-ray  aatronomy;  the  effecta  of  nuclear  axploalcMia, 

I magnetic  atorme,  and  aolar  activity  on  the  earth'a  atmoapharo,  ionoaphara,  and 

t magnetoaphera;  ^a  effecta  of  optical,  electromagnetic,  and  particulate  radio- 
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